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Abstract 

Controllability of chemical reactions is an important problem in 
chemical engineering science. In control theory, analysis of the con¬ 
trollability of linear systems is well-founded, however the dynamics of 
chemical reactions is usually nonlinear. Global controllability proper¬ 
ties of chemical reactions are analyzed here based on the Lie-algebra 
of the vector fields associated to elementary reactions. A chemical 
reaction is controllable almost everywhere if all the reaction rate co¬ 
efficients can be used as control inputs. The problem where one can 
not control all the reaction rate coefficients is also analyzed. The re¬ 
action steps whose reaction rate coefficients need to be control inputs 
are identified. A general definition of consecutive reactions is given, 
and it turns out that they are controllable almost everywhere by using 
the reaction rate coefficient of only one reaction step as control input. 

Keywords: kinetics, chemical reactions, controllability, Lie-algebra 
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1 Introduction 


Controlling chemical reactions is a key issue in modern chemical engineering 
science, see e.g. na, where control of nonlinear chemical reactors is consid¬ 
ered. The well founded methods in control engineering are based on linear 
dynamical systems, however the dynamics of chemical reactions is usually 
nonlinear. Control and controllability of nonlinear systems is only available 
for a relatively small class of systems, and is more involved than control of 
systems with linear dynamics, see e.g. dDlE]. Controllability is a fundamen¬ 
tal property of dynamical systems, since it tells us whether it is possible to 
reach the desired goal or not. 

Controllability of chemical reactions is usually analyzed using control the¬ 
ory developed for linear systems, and the linear model is acquired by lineariz¬ 
ing the dynamics at an operating point. Controllability of chemical reactions 
that have a positive equilibrium is established in [5j based on the linearized 
dynamics. In [2D], controllability analysis of a liquid-phase catalytic oxida¬ 
tion of toluene to benzoic acid was done based on linearization at hve different 
operating points, and controllability analysis of polymerization in different 
operating points was done in m Controllability analysis of protein glyco- 
sylation was done based on a linear model identihed from measurements in 
[T]. The connection of controllability and structure of chemical reaction was 
analyzed in [T3, while controllability of chemical processes was analyzed in 

mmmmm- 

Controllability analysis in operating points however only gives local re¬ 
sults. The application of Lie-algebra rank condition was suggested for the 
controllability analysis of chemical reactions in |3], however no general re¬ 
sults were given. A review of controllability analysis of chemical reactions 
was done in |21) . 

We analyze the Lie-algebra of the vector helds related to the reaction 
steps, and use the Lie-algebra rank condition as suggested in [1] to get global 
controllability results for chemical reactions. We consider the reaction rate 
coefficients of the reaction steps as control inputs and seek to hnd the lowest 
number of control inputs needed to control the system. The analysis done in 
this paper is symbolical, contrary to the analysis done in the literature that 
is usually numerical with some few exceptions, e.g. [5]. Our investigation 
will be restricted to those reactions for which the stoichiometric and kinetic 
subspaces coincide [7]; however this is not a too strong restriction. 

We give the model structure of kinetic equations that is used to model 
chemical reactions, and the Lie-algebra rank condition in Section [2j We 
show that chemical reactions are controllable almost everywhere if the con¬ 
trol inputs are the reaction rate coefficients of the reaction steps in Section 
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[21 In Section 01 we further analyze the Lie-algebra of the vector helds cor¬ 
responding to reaction steps, and give results on controllability when not all 
the reaction rate coefficients are control inputs. We give a lower and upper 
bound for the required number of control inputs, and identify the reaction 
steps whose reaction rate coefficients need to be control inputs. We dehne the 
initializer reaction steps and prove that their reaction rate coefficients need 
to be control inputs, that has already been shown by experiments in m on a 
polymerization example. We dehne the consecutive reactions, and show that 
they can be controlled with only one control input. We give two examples 
for consecutive reactions in Section |5l and analyze their controllability and 
region of controllability. The paper ends with a discussion in Section [6l 

2 Kinetics of reactions and controllability 

2.1 The differential equation of reaction kinetics 

Let the number of species participating in a chemical reaction be M, and the 
species be denoted by X^, X 2 ,..., Xj^. Suppose, that the reaction consists of 
R reaction steps. We call the linear combination of the species on the left- 
hand side of a reaction step the reactant complex, and the linear combination 
of the species on the right-hand side of a reaction step the product complex, 
so a reaction step is dehned by giving the reactant complex that transforms 
into a product complex with a given reaction rate coefficient. Let a{m,r) 
denote the number of species X^^ in the reactant complex of the rth reaction 
step and let /3{m, r) denote the number of species X^^ in the product complex 
of the rth reaction step. Let be the reaction rate coefficient of the rth 
reaction step. The chemical reaction system (or simply: reaction) is dehned 
by the R reaction steps |3] 

M M 

^ a(m, r)X^ r)X^, r = 1,2,...,R. (1) 

m=l m=l 

The matrix 7 = /3 — ck is called the stoichiometric matrix, and its columns 
are called the reaction step vectors. The usual mass action type differential 
equation of the reaction ([T]) is 

R 

x{t) = r)a;(f)“('’’’) (2) 

where Xm{t) is the concentration of the species X^^ at time t, x{t) is the 
vector of concentrations with its mth element being 7 (',t) denotes 
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( 3 ) 


the rth column of the matrix 7 , and is the monomial 




— d,2 


"f'M 


Definition 1 We call the image space of the stoichiometric matrix 7 the 
stoichiometric space. 

The stoichiometric space is the tangent space of the system ([2]) in almost all 
practically relevant cases. As to the exceptions, which we disregard here, see 
[7]. Since the stoichiometric space is the tangent space of the system ([2]), the 
dynamical system is only able to change its state locally in the sets of the 
form (a^o + Ran 7 ) fl called reaction simplexes [3]. 

2.2 Controllability of nonlinear systems and the Lie- 
algebra rank condition 

Consider an input affine system with the differential equation and initial 
condition 


x{t) = + '^gj{x{t))uj{t), a:;(0) = Xq e (4) 

i=i 

We say that the system is smooth if f and q,, On,..., o r are smooth func- 

Definition 2 System (jl]) is called controllable at a point x* G if there 
exists an open neighborhood U of x* such that if ®(0) G U then there exists 
Ui{t), U 2 it ),..., uj{t),t > 0 such that lim^^+oo x{t) = x*. 

Note that we only deal with controllability, and do not discuss the problem 
of finding the input functions that are required to get into the desired state. 

A conventional way of analyzing controllability of nonlinear smooth input 
affine systems like (jl]) is based on the Lie-algebra generated by the control 
vector fields Qj, j = 1, 2,..., J and the Lie-bracket of the drift vector field / 
and the control vector fields. 

Definition 3 The Lie-bracket of two smooth vector fields c^, t/? G M^) 

is 

[c^, -0] = (5) 
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Definition 4 A Lie-algebra jl] is a vector space V over the field F equipped 
with the binary operation [•, •] called the Lie-bracket (or Lie-product), if it 
satishes the following properties. 

1. Bilinearity, i.e. Vci, C 2 G F and V®, y, z E V 

[cix + c 2 y,z] = ci[x,z] +C 2 [y,z] (6) 

[z,cix + C 2 y] = Cl [ 2 , *] -f C 2 [z, y] (7) 

2. \/x G V, [x, x] = 0. 

3. The Jacobian identity, i.e. \/x, y, z E V 

[x, [y, ^]] + [y, *]] + y]] = O- (8) 

The hrst two properties of the Lie-bracket imply anticommutativity, i.e. 
'^x,y E V, [x,y] = —[y,x]. In this article, we will use the Lie-algebra on 
the vector space of smooth functions M^) over the held M equipped 

with the Lie-bracket dehned by ([S]). 

Definition 5 The Lie-algebra C generated by the vector helds and / 2 , 
denoted by £ = Lie{fi,f 2 ) is the smallest subspace of that 

satishes the following properties |2]: 

1- /i,/2e£; 

2. li x,y E C, then [a:;, y] E C. 

Definition 6 The object that assigns to every point x E a linear sub¬ 
space of is called a distribution nni. Distributions are usually denoted 
by the symbol A. 

Definition 7 A vector held / is said to be an element of a distribution A 
(f G A), if the value of f is in the subspace assigned to A in every point, 
i.e. Wx E fix) E A(x). 

Definition 8 A distribution Ac is called the controllability distribution ng 
of the system (jl]), if 

1. Wx E span{g^ix),g 2 ix),.. .,gjix)} C Ac. 

2. It is invariant under the vector held /, i.e. if r G Ac, then [r, /] G Ac- 

3. It is involutive, i.e. if ti,T 2 G Ac, then [ti,T 2 ] G Ac. 
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The above conditions mean that the controllability distribution is spanned 
by the vector helds from the Lie-algebra generated by the vector helds q, 
and j = 1 , 2 ,..., J |2]. 

Definition 9 The controllable subspace of the system (0]) is the linear combi¬ 
nation of all the space variable vectors that can be reached using appropriate 
inputs Ml, U 2 , ■ ■ ■, uj. 

Theorem 1 The dimension of the subspace Ac (a?*) is the dimension of the 
controllable subspace of the system at the point x*, so the system is control¬ 
lable in X* if and only if dimAc(a?*) = M. 

This is also called the Lie-algebra rank condition Note that for linear 

systems with the differential equation 

x{t) = Ax{t) + Bu{t) (9) 

the controllability distribution assigns the same subspace for every point, and 
this subspace is the image space of the controllability matrix 

Me =[ B AB A^B ... A^-^B ] . (10) 

So a linear system is controllable if and only if rank Afc = M. For example, 
in Mathematica the function ControllableModelQ can be used to determine 
if a linear system is controllable or not. 

3 Controllability of chemical reactions if all re¬ 
action rate coefficients are control inputs 

As a hrst step we consider chemical reactions in the form of ([2]) assuming 
that all of the reaction rate coefficients can be varied independently, so all 
reaction rate coefficients are control inputs. The chemical reaction under 
consideration thus is the driftless smooth input affine system 

R 

x{t) = 7 (-, r)x{t)"^'’^^Ur{t) (11) 

r=l 

that is obtained from ([2]) with the substitution kr = Ur{t) for r = 1,2,, R. 
This chemical reaction consists of R reaction steps with M species. Usually 
the M differential equations in (ITT]) are not linearly independent, in spite of 
the fact that usually R is much larger than M; in combustion chemistry, the 
Law’s law [12] says that R ~ 5M. 
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Since the trajectories of a chemical reaction system are in the reaction 
simplex that is contained in Ran 7 , the dimension of the controllable subspace 
can not be larger than the dimension of the image space of 7 . 

Definition 10 We say that a chemical reaction is controllable at a point 
X* G if dim Ac(x*) = rank 7 . 

The controllability of chemical reactions in the form of fllll) have already 
been analyzed in [5], where the author proved the following two theorems: 

Theorem 2 Suppose that the system ( 1111 ) has a positive equilibrium x* for 
Ui = U 2 = ■ ■ ■ = Ur = 1. Then for the control matrix B of the linear system 
(|9|) obtained from (ITT]) by linearization at the equilibrium point x*, we have 
that ranks = rank 7 , i.e. the linearized model is controllable, meaning that 
the system is controllable at the equilibrium point x*. 

Theorem 3 Assume that a chemical reaction has a positive equilibrium x* 
and that rank 7 = M. Then the system is controllable in the positive orthant. 

We prove that chemical reactions in the form of (fTTj) are controllable 
almost everywhere without the need to assume the existence of a positive 
equilibrium; moreover they are controllable in the positive orthant. 

Theorem 4 A chemical reaction with all reactions rates as control inputs 
is controllable almost everywhere. Moreover, it is controllable at every point 
of the positive orthant. 

Proof. The controllability distribution is spanned by the vector helds from 
the Lie-algebra generated by the vector fields gi,g 2 , ■ ■ ■, gR, and the vector 
fields [/,gi],[/,g 2 ],---so it is true, that 

span {^ 1 , c/ 2 ,..., c/r} C Ac ( 12 ) 

and since for every strictly positive (or negative) vector x this subspace is 

span {c/i, ^ 2 , ■■■,9r} = span { 7 (-, l)ci, 7 (-, 2 )c 2 ,..., 7(-, R)cr} (13) 

with Cr = being a nonzero constant for a fixed a; > 0 (or cc < 0) 

for r = 1 , 2 ,..., R, this implies that Ac equals to the image space of 7 , so 
dim Ac = rank 7 . Therefore, if none of the coordinates of x is zero, then 
the system is controllable. Note that it is possible that Xm = 0 for some 
m G { 1 , 2 ,..., M} and the system is still controllable. 

Another result from [5] is that if we make a reaction step reversible, and 
add this reaction step to the original reaction, then under some conditions, 
the controllability properties of the resulting system does not change. 


Consider a reaction that has a positive eqnilibrinm. Let the inpnt matrix 
of this reaction linearized in the positive eqnilibrinm be B. Take a reaction 
step, and make it reversible, and add this step to the original reaction. Snp- 
pose that the resnlting reaction also has a positive eqnilibrinm, and let the 
inpnt matrix of the linearized model be B. Then the rank of these matrices 
are the same |3], as it is stated in the following theorem. 

Theorem 5 Let ns assume that a given reaction has a positive equilibrium. 
If we make a reaction step reversible and the obtained reaction also has a 
positive equilibrium then rank.B = rank.B. 

We show, that if we make any reaction step reversible, then it does not change 
the dimension of the controllability distribution regardless of the existence 
of a positive equilibrium. 

Theorem 6 If we make any reaction step reversible, then the resulting 
chemical reaction is controllable almost everywhere; moreover, it is control¬ 
lable at every point of the positive orthant. 

Proof. Make the rth reaction step reversible, i.e. the reaction step vector 
of the added reaction step is ct{-,r) — f3{-,r) = —j(-,r) and the reactant 
complex of the new reaction step is defined by the vector /3(-,r), so the 
differential equation of the new reaction step is 

(14) 

The 7 stoichiometric matrix of the extended reaction is 

7=(7(-,l) 7 (-, 2 ) ... 7 (-,r) ... - 7 (-,r) ) (15) 

so it immediately follows that rank 7 = rank 7 . The Ac controllability 
distribution of the extended reaction possesses the property, that 

span { 7 (-, l)ci, 7 (-, 2 )c 2 ,..., 7(-, ^)cr, -7(-, r)cn+i} C Ac (16) 

where the left-hand side is the image space of 7 if the constants Ci, C 2 ,..., cji 
are not zero. 

Based on the results of this section, we can conclude that chemical reac¬ 
tions are controllable at almost every point if all the reaction rate coefficients 
can be varied independently to control the system. Moreover, controllabil¬ 
ity is not harmed even if we make any reaction step reversible. In the next 
section we analyze the practically more important case when not all reaction 
rate coefficients can be used as control inputs. 
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4 Controllability of chemical reactions with a 
reduced number of control parameters 


Suppose, that we can not control all the reaction rate coefficients indepen¬ 
dently, only a subset of them. Suppose without the loss of generality, that 
there are K number of controllable reaction rate coefficients, and the reaction 
steps are indexed so that the reaction rate coefficient of the first K reaction 
steps can be controlled. Let the index set of the controllable reaction rate 
coefficients be denoted by / = { 1 , 2 ,..., iL}, while the index set of the re¬ 
maining reaction rate coefficients be denoted by iV = {K -|-1,..., R}. In this 
case 1 < K < R. Let the vector field corresponding to the rth reaction step 
be denoted by Qr, i.e. 

c/^ : X (17) 

then the differential equation of the reaction with K control inputs is 


= '^kngn{.x{t)) + ^g^{x{t))Ui{t). 


n&N 




i&I 


(18) 


Definition 11 The controllability distribution of the reaction flTSl) with K 
control inputs is the distribution if 


1. Va) e span (a;), c/2 (a;),..., (a;)} C A[f^(a;). 

2. It is invariant under the vector field /, i.e. if r G A[^\ then [r, /] G 
a(^) 


3. It is involutive, i.e. 


if Ti, Ta G AJ^^ 




then [ti, t ^] G A[^^. 


The controllability distribution A[^^ is thus spanned by the vector fields 
from the Lie-algebra generated by the vector fields gi,g 2 ,- ■ ■ ,gK the 
vector fields [/,g'i], [f,g 2 ], ■ ■ ■, [/ifl'ic]- Therefore in the following analysis, 
we will need the Lie-brackets of the vector fields corresponding to reaction 
steps of the form of fITTj) . We will show, that at each fixed x G the Lie- 
bracket of the vector fields corresponding to two reaction steps is the linear 
combination of the reaction step vectors. First, we introduce the following 
notation 




0 , if a(m, r) = 0 , 

^a{-,r)-em otherwise. 


(19) 


with Em being the mth M-dimensional unit vector. Defining the diagonal 
matrix Dr{x) = diag {dl{x), ^ 2(^)5 • • • 5 can write the derivative 
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of the monomial in the short form 

^ a{l,r)d[{x) a{ 2 ,r)d 2 {x) ... a{M,r)d'^j^{x) ) 

= a{-,ryDr{x). (20) 

Utilizing these notations, we may formalize the Lie-bracket of vector helds 
corresponding to reaction steps as follows. 


Lemma 1 Let QijQj be vector fields in the form of flTTli . Then their Lie- 
bracket at the point x G is 

[9i,9j]ix) = 'y{-,j)oL{-,jyDj{x)j{-,i)x^^'’^^ ( 21 ) 

-7(-, i)Q:(-, i)^Di{x)j{-J)x‘^^'’^\ 


Proof. Since the Lie-bracket of two vector fields is = g'j9i — g'iQj by 

()S|) , and the vector helds are of the form of flTT)) , we only have to analyze the 
derivative of the vector helds. The derivative of the vector held g^ at a point 
X G is 


/ digi^i{x) d2gi,i{x) ... dMgi,i{x) \ 

digi^2{x) d2gi,2{,x) ... dMgi,2{x) 

9^{x) = . . 

\ digi^M^x) d2gi,M{x) ... dMgi,M{x) 

where Oimix) = 7(m, i)a3“('’*d 1 


( ( 7 ( 1 , ^ 

( 7 ( 2 , z)a;“f 

V (7(M,i)a;“('’*))' / 
1221 




9\{x) = 7(-,i)Q:(-,i)^A(*), 


(23) 


where we have utilized the notation introduced before the lemma. Since the 
vector held at a point x G is gj{x) = 'y{-, j)x°'^'’^\ the statement of 
the Lemma follows after substitution into the dehning equation ([5]) of the 
Lie-bracket. 


Lemma 2 Wx G Di{x)j{-, j) = 0 77 a{m,i)''y{m, j) = 0, m = 

1,2,...,M. 

Proof. Expanding the term Di in the following expression, we get 

a(-,i)^A(®)7 (-,j) = a{l,i)-f{l,j)d\{x) + a{2,i)-^{2J)d\{x) + 

+ ... + a(M,i)7(M,j)d)v^(a;) (24) 

where df^fx) is either the zero polynomial, or the derivative of a hxed poly¬ 
nomial with regard to its mth variable, for any m. However A = 0 for all 
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m G {1, 2,..., M} if and only if q:(-, i) = 0, in this case the statement holds 
trivially. So snppose that q:(-, i) ^ 0. If = 0 for some m, then a{m, i) = 0 
by the dehnition of so a(m, i) = 0 for that m. So we can snppose 

withont the loss of generality that d^ 7 ^ 0 for any m G {1, 2,..., M}. Since 
the nonzero polynomials d^, m = 1,2,..., M are all the derivatives of the 
monomial it follows that they are linearly independent, i.e. 

M 

a{m, j)dl^{x) = 0^ a{m, j) = 0, m = 1,2,..., M (25) 

m=l 

that conclndes the proof. 

If we use the notation © for the elementwise (or Hadamard) product of 
two vectors, the condition for Lemma[2]can be written as q:(-, i)Q^{-,j) = 0. 

Definition 12 Let ©j : x 0 !.(-,i)^Di{x)j(-, j). 

Using this notation, the Lie-bracket of two vector helds of the form flTTj) is 

[9i,9j]ix) = Kj,i{x)x^^'’^'^-r{-,j) - Kij{x)x^^'’^^j{-,i). (26) 

In the analysis of the Lie-bracket of vector fields corresponding to reaction 
steps, species that participate in both the reactant and product complex of 
a reaction step by the same amount do not count. We call this species the 
direct catalyst of a reaction step. 

Definition 13 The mth species is the direct catalyst of the rth reaction 
step, if the following two (equivalent) statements hold. 

1 . 7 (m, r) = 0 , but a{m,r) 7 ^ 0 . 

2. The mth species participates in the reactant and product complex of 
the rth reaction step by the same amount. 

Note that in the reaction step 

X + 2Y ^ 2X +Y 

neither the species X nor the species Y is a direct catalyst according to our 
definition. 

Definition 14 The mth species is the direct catalyst of a reaction, if the 
following two (equivalent) statements hold. 

1. 7(m, •) = 0. 
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2. If the mth species participates in a reaction step, then it is a direct 
catalyst of that reaction step. 

Lemma 3 The Lie-bracket of two vector helds corresponding to reac¬ 
tion steps written in the form of ffT7|) at the point x G is 

1. [^j, gj\ {x) = 0, if there exists no species in the reaction that participates 
in both the zth and jth reaction step or if it participates, it is a direct 
catalyst of the reaction steps. 

2. [gi,gj]{x) = —Kij{x)x"^'’^'>j{-,i), if there exists at least one species 
that participates in the reactant complex of the ith reaction step, and 
participates in the jth reaction step, but not as a direct catalyst of the 
jth reaction step, and there exists no species that participates in the 
reactant complex of the jth reaction step but does not participate in 
the zth reaction step except as a direct catalyst. 

3. [gj^,gj](x) = +Kj^i{x)x°'^'’'^^j(-, j), if there exists at least one species 
that participates in the reactant complex of the jth reaction step, and 
participates in the fth reaction step, but not as a direct catalyst of the 
fth reaction step, and there exists no species that participates in the 
reactant complex of the ith reaction step but does not participate in 
the jth reaction step except as a direct catalyst. 

4- [gi,gj]ix) = Kj^i{x)x°‘^'’^^j{-,j) — Kij{x)x°‘^'’^'>^{-,i), if there exists a 
species that participates in the reactant complex of the ith reaction step 
and participates in the jth reaction step but not as a direct catalyst, 
and there exists a (possibly different) species that participates in the 
reactant complex of the jth reaction step and participates in the zth 
reaction step, but not as a direct catalyst. 

Proof. (1): The condition is equivalent to cx{-,i) = 0 and Q:(-, j) O 

7 (-,f) = 0, which implies that [gi,gj] = 0 because of Lemma [2j 

(2) : The statement is equivalent to cx{-,i) © 7(-, j) ^ 0 but 0 

7 (-,z) = 0, so by Lemma[2]Kjj ^ 0 but Kj^i = 0, that implies [gi,gj]{x) = 
—Kij{x)x°‘^'’^^^{-,i) because of fl26|) . 

(3) : The consequence of (2) and the anticommutativity of the Lie-bracket. 

(4) : The consequence of (2) and (3) and the bilinearity of the Lie-bracket. 

In order to prove our next theorem, we need the following lemma: 

Lemma 4 Let P{x) = ( f^{x) f^ix) ... f^{x) ) be an n X n matrix 
whose elements are polynomials of cc G M” of hnite order. If there exists 
a point X* G M"' such that det P(a^*) ^ 0, then det P ^ 0 holds almost 
everywhere in M"'. 
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Proof. Since the elements of the matrix are hnite order polynomials, the 
determinant of P is also a hnite order polynomial. Since there is a point at 
which det P 7 ^ 0 , this polynomial is not the zero polynomial, thus the set of 
the roots of the determinant has zero measure in M", so detP 7 ^ 0 almost 
everywhere in M"'. 


Theorem 7 If a chemical reaction is controllable at a point x* G with 
the inputs K, then it is controllable almost everywhere. 

Proof. If the reaction is controllable at a;*, then dirnA^^^ = rank 7 , which 
means that there are rank 7 number of independent vectors spanning 
The controllability distribution a[^^ is spanned by the vector helds from 
the Lie-algebra generated by vector helds of the form of flT7|) . which are 
polynomials with each monomial being multiplied with one of the reaction 
step vectors because of Lemmas [I]-[21 thus by collecting the reaction step 
vectors we get that 

span{pi(a;)7(-,l),p2(a;)7(-,2),...,pK(a;)7(-,^)} ^ (27) 


with pr, r = l,2,...,P being some polynomials. If the reaction is control¬ 
lable at a point x*, then 

R 

'^CrPr{x*)-f{-,r) = 0 (28) 

r=l 


if and only if there are rank 7 number of c^s that are zero. Let rank 7 = S. 
Suppose without the loss of generality, that the hrst S component of the vec¬ 
tors pi{x*)j{-,l), p2{x*)j{-,2),..., ps{x*)j{-, S) are linearly independent. 
Denote the choice of the hrst S components of the polynomial vector Pi^{-, i) 
by Pi^i', Then linear independence can be expressed so that the system 
of equations dehned by 


( pi7(-,l)[®](a;) p2j{-,2y^^{x) 


Psl{-,Sy^^{x) ) 


/ Cl \ 

C2 

\cs J 


= 0 (29) 


at the point x* has only the trivial solution. This yields that the determinant 
of the matrix on the left-hand side of (|29|) is not zero, and by Lemma IU the 
determinant of that matrix is not zero almost everywhere in M'^, so the vectors 
are linearly independent almost everywhere. 

In the following proof we give an upper estimate for the required number 
of control inputs. 
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Theorem 8 A reaction is controllable almost everywhere with rank 7 num¬ 
ber of reaction rate coefficients used as control inputs, i.e. K = rank 7 . 

Proof. Choose as control inputs the reaction rate coefficients of reaction steps 
whose reaction step vectors are linearly independent. Since the controllability 
distribution contains these vector fields corresponding to rank 7 number of 
independent reaction steps, the dimension of the distribution is rank 7 almost 
everywhere. 

In practice, usually rank 7 < R, so controlling only rank 7 reaction rate 
coefficients may be much easier than controlling the reaction rate coefficients 
of every reaction. The question is, whether we can further decrease the 
number of control inputs. 

Definition 15 Let the reaction rate coefficients of the hrst K reaction steps 
be control inputs. A reaction step is said to be the critical reaction step 
of the control system with the K inputs, if removing the rate coefficient 
of that reaction step from the control inputs reduces the dimension of the 
controllability distribution almost everywhere. 

Simply saying, if we remove the reaction rate coefficient of a critical reac¬ 
tion step from a controllable system then the resulting system will not be 
controllable. Identihcation of critical reaction steps is crucial in hnding the 
required control inputs. 

Definition 16 A reaction step is said to be an initializer, if the species of its 
reactant complex do not participate in any other reaction steps, except if such 
a species is a direct catalyst in the given step, and the species participating 
in the reactant complex are not all direct catalysts. 

For example, in the reaction 


X + U X + V (30) 

X + Y-^X + Z (31) 

both reaction steps are initializers, while in the reaction 

X + U U +V 

X +V x +w 

only the first reaction step is an initializer. 

Initializers have key role in the control of chemical reactions, because they 
are always critical reaction steps. 


(32) 

(33) 
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Theorem 9 An initializer is a critical reaction step for all K input sets that 
contain the initializer. 

Proof. Let the ith reaction step be an initializer. Due to the definition of 
the initializer, there exist an m G { 1 , 2 ,... ,M} such that ^ 0 , but 

7 (m, j) = 0 for all j G {1, 2,..., R}, j ^ i, so the reaction step vector corre¬ 
sponding to an initializer is linearly independent from all the other reaction 
step vectors of the chemical reaction. 

Suppose indirectly, that there exists a number K such that has 

dimension rank 7 almost everywhere, but i ^ I = (1,2,...,iL}. We will 
show, that this can not happen. 

Since i ^ I, and is linearly independent from other column vectors 

of 7 , 7 (-,f) is linearly independent from all k ^ I almost everywhere. Let 
Q denote the Lie-algebra generated by the vector fields Qf., k E I. Due to 
Lemma [3l this Lie-algebra consists of vector fields that can be written as 

g = ^li-,k)pk (34) 

k£l 

with pk being some polynomials. At almost every x G the vector g{x) 
is linearly independent of since it is the linear combination of vectors 

7 (-, k),k E I that are linearly independent from 

Let us examine the vector fields of the form [f,G]. Because the Lie- 
algebra is bilinear, this can be written as 


[f,S] 


[kiQi, Q] + 


X] knQn, G 
jiGN^n^i 


(35) 


Vector fields of the form knOn, G] are linearly independent of 7 (-, i) 

almost everywhere due to the same argument as before, since the vectors 
E N,n ^ i are linearly independent from almost everywhere. 

Albeit gi{x) = 7 (-,is linearly dependent of 7 (-,z), the vector 
fields [kig^^Q] are also linearly independent of since the Ah reaction 

step is an initializer, i.e. there exists no species that is in the reactant complex 
of the Ah reaction step, but participates in the other reaction steps except as 
a direct catalyst, so by Lemma [3l the Lie-bracket of these vector fields does 
not contain a component linearly dependent of 7 (-,f). 

The Lie-algebra generated by Q and [/, Q] is also linearly independent 
from 7 (-,f) almost everywhere, since the vector fields in the generators are 
linearly independent from almost everywhere. 

The controllability distribution is formed by vector fields from the Lie- 
algebra generated by Q and [f,G], but 7 (-,f) is linearly independent from 
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iK^ 

these vector fields, thus dim ^ < rank 7 almost everywhere, which implies 
the system is not controllable, that is a contradiction. 

This theorem states that the reaction rate coefficient of an initializer needs 
to be controlled in order to control the chemical reaction, so a lower estimate 
of the required number of control inputs is the number of initializers in a 
reaction. Note that the number of initializers can be easily found utilizing 
the ex and 7 matrices. Usually the required number of control inputs is higher 
than the number of initializers, but there are some special reactions that can 
be controlled by controlling the reaction rate coefficients of its initializers. 
The following are also true about initializers: 

1. A step of a reversible pair can not be an initializer. 

2. If the reaction is weakly reversible, then there is no initializer, but it is 
possible to find initializer classes, which we define later. 

3. If any of the species of the reactant complex are zero initially, then the 
initial step can not start, and the reaction is not controllable at the 
initial point. 

Definition 17 A reaction is called a consecutive reaction if it has the fol¬ 
lowing properties. 

1. It has exactly one initializer. 

2. If we remove the initializer, the resulting reaction has exactly one ini¬ 
tializer (or it is empty), so it has the previous property. 

3. If we remove the initializer, the resulting reaction has the previous two 
properties. 

In consecutive reactions the reaction steps follow each other serially, so they 
may also be called serial reactions. The examples in Section [5] will show that 
our definition of consecutive reactions is more general than the one usually 
used in chemical kinetics. A consecutive reaction has only one initializer, 
and the reaction is controllable by controlling the reaction rate coefficient of 
the initializer only. 

Theorem 10 A consecutive reaction is controllable with only one input, if 
and only if the input is the reaction rate coefficient of its initializer. 

Proof. If the reaction is controllable with one input, then the input must be 
the reaction rate coefficient of the initializer due to Theorem O 
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Now we prove, that if the control input is the reaction rate coefficient 
of the initializer, then the reaction is controllable. Let the initializer of 
the consecutive reaction have the index ii. If we remove the initializer, we 
get a consecutive reaction due to the dehnition of consecutive reactions, or 
an empty reaction. Let the initializer of the new reaction be the reaction 
step that has the index *2 in fhe original reaction. The reactant complex of 
the f 2 th reaction step has at least one species that participates in the iith 
reaction step and is not a direct catalyst in either reactions, otherwise the 
original reaction would have the i 2 th reaction step as an initializer as well. 
The reaction step with index ii does not have common (non direct catalyst) 
species with reaction steps other than the one with index i 2 , since if it had a 
common non direct catalyst species with an reaction step with index j ^ i 2 , 
it would mean that removing ii would result in a reaction with both ^2 and 
j being initializers (or none of them), that contradicts the dehnition of a 
consecutive reaction. So the reaction step with index ii has common species 
only with reaction step with index * 2 - Remove the reaction steps with indices 
ii and * 2 - The resulting reaction is either empty, or it has one initializer, let 
this be the reaction step with index in the original reaction. Using similar 
argument as before, the reaction step with index *2 has a non direct catalyst 
species that participates in the reactant complex of the reaction step with 
and the reaction step with index *2 does has common species with reaction 
steps with indices R and only (except for species that are direct catalysts 
of the f 2 th reaction step). We hnd the i„th reaction step by removing the 
hrst n — 1 initializers recursively. Then the fnfh reaction step has common 
species (other than direct catalysts) only with reaction steps with indices 
in-i and in+i, and the reactant complex of the reaction step with index in 
has a common non direct catalyst species with the reaction step with index 
1 • 

The indices of the reaction steps of the consecutive reaction can be ar¬ 
ranged into the set {ii,i 2 ,, in}, and the drift and the control vector helds 
of the consecutive reaction can be written as 

R 

f : 

r=2 

g ■ (36) 

The controllability distribution consists of the vector helds from the Lie- 
algebra generated by the vector helds g and [f,g]. The vector held g spans 
7(-,fi) almost everywhere because of its dehnition. The Lie-bracket [f,g] 
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can be written as 


[f,9] = [di^Qn] + 


J^2 


(37) 


where the second term is zero due to Lemma O since the reaction step with 
index ii has common non direct catalyst species only with the reaction step 
with index 12 - The reaction step with index 12 has a non direct catalyst 
species in its reactant complex that participates in the iith reaction step not 
as a direct catalyst, but there is no species in the reactant complex of the 
iith reaction step that participates in the i 2 fh reaction step as a non direct 
catalyst, so the application of Lemma [3] results in the vector held at 

the point x G written as 


[f.9\{x) = {x) = (38) 

so g and [f^g] span { 7 (-, ii), 7 (-, * 2 )} almost everywhere. Due to the dehni- 
tion of the vector held [f^g] can be interpreted as a vector held corre¬ 
sponding to linear combination of hctive reaction steps with the same reaction 
step vector as 7 (-,Z 2 ), and the reactant complex vectors being element-wise 
less than or equal to a.{-, * 2 ) + Q:(-, * 1 ), so the species in the reactant complex 
vector of the hctive reaction steps are from the species of the reactant com¬ 
plexes of the reaction steps with indices 7 and 12 , that does not participate 
in any other reaction steps. So the hctive reaction steps related to [f,g] 
has common species with the reactant complex of the isth reaction step, and 
the reactant complex of the hctive reaction step related to [fyO] has common 
species with the iith and i 2 fh reaction step (except for direct catalysts). This 
implies, that the vector held [/, [/,fl']] can be written as 


[/, [/, 9 ]] = [9i2, «i2,n7(-, * 2 )] + [9i3, «i2,ii7(-, * 2 )] (39) 

that can be written at a point x G as 

[/, [/, 9 ]]{x) = pi, 2 (a;) 7 (-,^ 2 ) + P 2 , 2 (a;) 7 (-,d) (40) 


where pi _2 and ^ 2,2 are polynomials in which the monomials divide either 
Q]. vector held [/, [f,9]] can be interpreted as a vector 

held corresponding to the linear combination of hctive reaction steps with 
reaction step vectors 7 (-, 12 ) and 7 (-, i^) and the species in the reactant com¬ 
plex of the hctive reaction steps are the species from the reactant complexes 
of the reaction steps with indices ii,i 2 ,i 3 . So the hctive reactions related to 
the vector held [/, [f,9] does not have a species in their reactant complex 
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that participates in the reaction steps with index v,r > 3 as a non direct 
catalyst, but there is a species in the reactant complex of the Z 4 th reaction 
step that participates in the Active reactions related to [/, [/,g']]. 

Continuing the same argument, we get that application of the Lie-bracket 
with the vector field f to the vector field g n — 1 times results in the vector 
field with the value 


n—1 

in) ( 4 ^) 

in the point x G This shows, that application of the Lie-bracket with 
vector field / to the vector field g n — 1 times results in a vector field that 
spans {7(-, * 2 ), 7(-, * 3 ), • • •, 7(', in)} almost everywhere if n < i?. This implies 
that the Lie-algebra generated by the vector fields g and [f,g] (and thus the 
controllability distribution) spans the subspace {7(-, 7 ), 7(-, * 2 ), • • •, 7(‘, 
almost everywhere, so it spans the image space of 7 almost everywhere, thus 
the dimension of the controllability distribution is rank 7 almost everywhere. 

It is possible, that a chemical reaction does not have an initializer. How¬ 
ever, in this case there is a group of reaction steps such that if they were 
united into one reaction step, the result would be an initializer. We will call 
these groups the initializer classes. 

Definition 18 We say that a group of reaction steps form an initializer class, 
if the species from the reactant complex of the reaction steps in the initializer 
class does not participate in the reaction steps that are not in the initializer 
class, except as direct catalysts. 

Theorem 11 There is a critical reaction step in each initializer class. 

Proof. The proof is indirect. Suppose, that the reaction rate coefficients 
corresponding to the reaction steps in the initializer class are not among 
the control inputs, and that the chemical reaction is controllable almost 
everywhere. Consider the Active reaction that is the linear combination of 
the reaction steps in the initializer class. This Active reaction can be used in 
the analysis, since the drift vector Aelds of the chemical system in question 
consists of the linear combination of the vector Aelds corresponding to the 
reaction steps on the initializer class and possibly other reaction steps and 
the Lie-bracket is bilinear. The linear combination of the reaction steps of 
the initializer class is an initializer by deAnition, that is a critical reaction 
step due to Theorem [9l so the chemical reaction is not controllable, that is 
a contradiction. 
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5 Examples 


Consider the following consecutive reaction with three species (M = 3) and 
two reaction steps {R = 2) 


X, ^ X 2 (42) 

2X2^X3 (43) 


This reaction has one initializer, it is the hrst reaction step with reaction rate 
coefficient ki. The matrices a,l3,^ formed by the reactant complex vectors, 
product complex vectors and reaction step vectors respectively are 



This is a consecutive reaction, so by Theorem[Tn]it is controllable with a single 
input if the control input is the reaction rate coefficient of the initializer, i.e. 
it is the reaction rate coefficient ki. The control and drift vector helds are 


g .x^ 


-1 

1 

0 


Xi, f ■. x\-^ k2 



(45) 


and the differential equation of the chemical reaction is 


x{t) 


f{x{t)) + g{x{t))u{t) 





(46) 


where x = (xi, X 2 , X 3 )''' and Xi,X 2 and X 3 is the concentration of the species 
X 3 ,X 2 and X 3 , respectively. The controllability distribution is spanned by 
the vector helds from the Lie-algebra generated by g and [f,g]- The vector 
held [/, g] at the point a; G is 




(» 7 )(®) - {f9){x) 

/-100W0\ /O 00 

10 0-2 k2xl - A;2 0 -4X2 0 

\000/\l/ \0 2x2 0 


( ° 

- 2 / 02 X 1 X 2 I -2 

V 1 


-2/02XiX27(-, 2). 



Xi 


(47) 
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The vector fields g and [f,g] are linearly independent almost everywhere, 
so the dimension of the controllability distribntion is at least two almost 
everywhere. Since the dimension of the stoichiometric space is rank 7 = 2, 
the reaction is controllable almost everywhere. 

Now we determine the set of nnll measure on which the system is not 
controllable. The system is trivially not controllable on the line Xi = 0, since 
in this case g{x) = 0. Now suppose, that Xi 7 ^ 0. The vector field [f,g] is 
zero at a ;2 = 0, however there may be other vector fields in the Lie-algebra 
generated by the vector fields g and [/, g] that are not zero at X 2 = 0 . 

The vector field [g, [fig]] at a point £C G is 


9. l/.all(a:) = (|/,9l'9 -9'l/,9l)(a:) = Ihxiixn -i,) 



(48) 


that is not zero, if 0:2 = 0, however it is zero, if xi = X 2 . So if xi 7 ^ 0 and 
a ;2 7 ^ 0 , then the stoichiometric space is spanned by the vector fields g and 
[f,9l and if Xi 7 ^ 0 and X 2 = 0, then the stoichiometric space is spanned by 
the vector fields g and [g, [/, g]]. This implies that the reaction is controllable 
everywhere except on the line xi = 0 . 

Consider the following consecutive reaction that consists of three reaction 
steps with four species: 


X, + X, Xj + X 3 

Xa^X^ 

-V" V" 

A4 —t A2 


(49) 

(50) 

(51) 


This reaction has one initializer, that is the reaction step with the reaction 
rate coefficient ki. The matrices cx.,13,^ formed by the reactant complex 
vectors, product complex vectors and reaction step vectors respectively are 



/I 

0 

0 \ 


( ° 

0 

0 \ 


/-I 

0 

0 \ 


1 

0 

0 

, /3 = 

1 

0 

1 


0 

0 

1 

a = 

0 

1 

0 

1 

0 

0 

, 7 = 

1 

-1 

0 


VO 

0 

1 


VO 

1 

oy^ 


V 0 

1 

- 1 / 


Let the control input be the ki reaction rate coefficient of the first reaction, 
and let the concentrations of the species X^, X 2 , X 3 and X 4 be denoted by 
Xh X 2 , x^ and 0:4 respectively. Then the control and drift vector fields are 





( M 


( 


0 


0 


1 

g : 

1 

X1X21 f : X k2 

-1 

X3 + ks 

0 


V Oy^ 


V W 




( 53 ) 
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and the differential equation of the chemical reaction system is 


f{x{t)) + g{x{t))u{t) 



( ° ^ 



k3X4.(t) 

1 

0 

-k2X3{t) 


1 

1 

CO 

B 


V 0^ 


Xi{t)x 2 {t)u{t). (54) 


We calculate the vector helds [f,g] and [/, [f,g]]. The vector held [f,g] 

is 


at the point a; G 

[f,9]ix) 


{ 9 '. 


)- 

-(/'»)( 

x) 



/ 

-X2 

- 

-xi 0 

0 \ 


( 

0 

\ 


0 


0 0 

0 


k3X4 


X2 


xi 0 

0 


-k2X3 

V 

0 


0 0 

^ J 


V 

k2X3 - 

ksXA / 


/ 0 

0 

0 

0 

\ 





0 

0 

0 

h 


0 



0 

0 

-k2 

0 


1 

XIX 2 



0 

k2 

-k3 

} 


V 



( 

— 

ksXA 

\ 




= a;i 


V 


0 

k2X2 + k'iXi 
-k2X2 


The vector held [/,[/, S']] at the point a; G is 

[/,[/,«/]](a?) = ([/,«/]'/-/'[/,c/])(a^) 

/ -k 2 kzXz + ^3X4 

k2k3X2 

— 

klx2 - k^Xi + k2k3X3 + 2k2k3X4 
\ -k2{k2X2 + k3X2 + 2k3Xi) 


(55) 


(56) 


The vector helds g, [f)g] and [f,[f,g]] are linearly independent almost 
everywhere, and they are in the Lie-algebra generated by the vector helds 
g and [f,g], so the controllability distribution contains these vector helds, 
and thus the dimension of the controllability distribution is three almost 
everywhere. Since rank 7 = 3, the chemical reaction system is controllable 
almost everywhere. 

Now we determine the set of points in which the system is not controllable. 
Since g{x) = 0 if xi = 0 or 0:2 = 0, it follows that the system is not 
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controllable in the plane Xi = X 2 = 0. Suppose now, that Xi 7 ^ 0 and 
a ;2 7 ^ 0. It can be easily verihed that the vector fields g, [f,g] and [/, [f,g]] 
are linearly independent in this case. So the system is controllable in every 
a; G except if a:i = 0 or X 2 = 0 . 

6 Discussion 

We have shown that if all the reaction rate coefficients of a chemical reaction 
are control inputs, then the chemical reaction is controllable almost every¬ 
where, and it is controllable everywhere in the positive orthant. We have 
shown, that making a reaction step reversible and adding it to the original 
reaction does not affect controllability. 

In chemical engineering practice, reducing the number of control inputs 
is of great interest. We have defined critical reaction steps of input sets, that 
must remain in the input set to keep controllability. However, the state of 
being a critical reaction step depends on the current input set. We have 
defined initializer reaction steps, and showed that they are always critical 
reaction steps independent of the input set. Thus the minimal number of 
control inputs can not be less than the number of initializer reaction steps. 
We have defined consecutive reactions, and proved that they can be controlled 
with a single input. 

It is possible, that a reaction does not have an initializer reaction step, 
however it is still possible to find an initializer class of the reaction steps. We 
have proved, that there is always at least one critical reaction step in every 
initializer class. So the minimal number of required control inputs is not 
less than the sum of the number of initializer reaction steps and initializer 
classes. 

Most of the results are specific for chemical reactions, and the theorems 
and definitions are expressed using chemical concepts like e.g. the species 
participating in the reactant and product complexes. However, Theorem [7] 
only uses the fact that chemical reactions are polynomial systems, so this 
theorem applies for any system whose dynamics is described by a differential 
equation with its right-hand side being a polynomial of the state variables, 
and being linear in the control inputs. So in the case of input affine polyno¬ 
mial systems, local controllability implies global controllability. 

It is a nontrivial question how the inputs are to be chosen in order to 
achieve a given final state, we shall return to this — important from the 
practical point of view — question later. 

In the future, we plan to investigate the observability of chemical reactions 
using the same approach. By using the Lie-algebra of the smooth vector 
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fields, we hope to find global results for observability to generalize the local 
results that can be found e.g. in [HI EJ- 
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